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Abstract. We carried out a NIR variable star survey toward the Large and Small Magellanic Clouds using the InfraRed Survey 
Facility (IRSF) at Sutherland, South African Astronomical Observatory. This survey project was initiated in Dember 2000, 
and since then, we kept monitoring a total area of 3 square degrees along the LMC bar and also an area of 1 square degree 
around the center of the SMC, sufficiently large to do statistical analysis and to make complete catalog of variable red giants 
in the Magellanic Clouds. The detection limits (S/N=10) of the survey are 17.0, 16.5 and 15.5 at J,H and K s , respectively. In 
this article, we present some results on infrared variables that are not detected by the previous optical surveys. We show that 
they do not fall on the standard period— K s magnitude relation for Mira-type variables pulsating in the fundamental mode. 
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INTRODUCTION 

Gravitational microlensing search projects (OGLE, ej* 
IToll ; MACHO, e.g., H; EROS, e.g., [j; MOA, e.g., Jl) 
observed the Magellanic Clouds and found a large num- 
ber of variable stars. These surveys were operated in the 
optical wavebands, and, most, if not all, of high mass- 
losing variables might have been missed by these sur- 
veys, because such stars are not readily detectable in the 
optical wave bands due to their surrounding circumstel- 
lar dust shells. 

To know pulsation properties of such infrared vari- 
ables, we carried out a near-infrared monitoring survey 
toward the Large and Small Magellanic Clouds using 
the IRSF at Sutherland, South Africa. The IRSF consists 
of a 1.4 m telescope and a near-infrared (NIR) camera, 
SIRIUS. SIRIUS can observe the sky in the three wave 
bands, J,H and K s simultaneously and each image has a 
field of view of about 7.6' x 7.6'. An area of three square 
degrees in the LMC and another an area of 1 square de- 
gree in SMC have been monitored since 2000. These ob- 
served areas are graphically illustrated in 16J] and H. The 
main goal of this survey is to detect infrared, very red 
variables that cannot be detected by optical surveys. 

OBSERVING STATUS 

Observing status can be summarized as follows. 

LMC Monitoring was started in Dec. 2000, and were 
continued until April 2008. Therefore, the baseline 
of the observing span is over 7 years (> 2500 days). 
About 80-90 observations were done during that 
period for the 3 square degrees monitoring area. 



SMC Monitoring was started in Dec. 2001, and were 
continued until October 2007. Therefore, the base- 
line of the observing span is over 6 years (> 2100 
days). About 100-1 10 observations were done dur- 
ing that period for the 1 square degree monitoring 
area. 



FINDING LIMITS AND OUR 
DISCOVERY SPACES 

The monitoring observations are made with 5 seconds 
x 10 dithering mode, which yield time-series imaging 
data of S/N = 10 at 15.5 mag in K s band. Then, how 
bright/faint variables did we detect in the survey? The 
bright limit is determined by the saturation, which is 
about 8 mag at K s band. The faint limit depends on the 
"pulsation amplitudes" of faint variable stars. We de- 
tected classical Cepheids with pulsation amplitude (AK S ) 
of about 0.3 mag (i.e., about 3 ff) at around 15 mag at 
K s band (corresponding to the pulsation period of a few 
days). 

Then, with these finding limits, what kind of variable 
stars can our infrared monitoring survey detect? Let us 
take the optical OGLE-II survey in a comparison with 
our infrared survey. In figure Q] we indicate our discov- 
ery spaces on the (I — K s ) v.s. Ks color-magnitude dia- 
gram (hatched areas). The horizontal lines indicate the 
saturation and detection limits of our survey. Therefore, 
any variable stars with AK S > 0.3 mag that locate be- 
tween these horizontal lines are expected to be detected 
by our survey, as well as the OGLE-II survey. The diag- 
onal lines are determined by combinations of the satura- 
tion and 10 C detection limits (about 13 and 19.5 mag 
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FIGURE 1. The (/ — K s ) v.s. K s color-magnitude diagram 
of variable stars in the LMC, showing our discovery spaces 
(hatched areas). Reverent data are taken from 18]. The arrow 
shows interstellar extinction vector for A v = 10 mag, using 
Weingartner & Draine (| 13]) model. 

at / band, respectively flTll ) of the OGLE-II survey and 
the 10 (7 detection limit of our survey. The segregation 
of discovery spaces between optical and infrared surveys 
are clear. Cepheid variables and moderately red variables 
are mutually detected by OGLE-II and our survey. Mean- 
while, only our survey can detect very red variables, that 
are fainter than 19.5 mag at /, but brighter than 15.5 mag 
at K s , and also bright variables that are brighter than 13 
mag at /, but fainter than 8 mag in K s (e.g., supergiants). 

Unfortunately, our survey could not detect extremely 
red objects, which are recently discovered by Spitzer 
space telescope in the LMC ([5]). Such extreme stars 
are embedded in their surrounding dust shells with quite 
large optical depths (the most extreme star is buried in 
the shell with T\^ m > 200!!). They are fainter than 16.3 
mag even at Ks band, which is fainter than our detection 
limit. 



TENTATIVE RESULTS 

Due to the long baseline of the observation, we detected 
long period variables, whose pulsation periods exceed 
a thousand days. Also, we found a number of infrared 
variables, which are not detected in the OGLE-II survey. 
Most of the infrared variables pulsate very regularly, 
and their light curves have sinusoidal-curve like shapes. 
Also, their pulsation amplitudes are large even at K s 
band (Ita et al. in preparation). Judging from these facts, 
it is likely that these infrared variables are Mira-type 
variables. 

It is well known that the period— luminosity relation is 
best defined at K band, where the atmospheric bands do 
not contribute much and the fluxes correspond to the pho- 
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FIGURE 2. Period— K s magnitude diagram of variable stars 
in the LMC. Small dots show OGLE-II optical variables, and 
big diamonds show some of the infrared variables (not all), 
that are not detected in the OGLE-II survey. The diagonal line 
show the period— K s magnitude relation for Mira-type variables 
pulsating in the fundamental mode, which is defined in Q]. 

tosphere. Ita et al. (J3l) analyzed OGLE-II light curves 
and obtained pulsation periods and amplitudes for more 
than 8000 optical variables. Together with their original 
K s band photometry, they made Period—^ diagram as 
shown in figure [2] Optical variables are shown in small 
dots. There are many sequences, and each of them corre- 
sponds to the different type of variable stars and also to 
different pulsation mode. See, for example, Wood ([ 14]) 
and Ita et al. ([7]) for explanations on the sequences. The 
big diamonds show some of the infrared variables de- 
tected in our survey. 

It seems that the K s band fluxes of the infrared vari- 
ables do not represent their photospheres any more. 
Wood (HI) studied MSX (0]) sources in the LMC, and 
found that very red Mira-like variables do not fall on 
the standard period— K magnitude relation for Mira-type 
variables (e.g., sequence C of Ita et al. [7], and it is shown 
by the diagonal line in figure |2). Infrared variables fall 
below the relation. Note that some of the infrared vari- 
ables seem to be on so-called Ag se quence D ( II 1410 Ah. 
However, it is just a matter of appearance. As mentioned 
above, the light curves of infrared variables are beauti- 
ful sinusoidal-curve like, and their AK S are large, which 
are very much differ from those of variable stars on "se- 
quence D". 

Here we assume that the reason of reduced K s band 
fluxes of the infrared variables are due to the circum- 
stellar extinctions from their surrounding dust shells. We 
also assume that all of the infrared variables are Mira- 
type variables pulsating in the fundamental mode, and 
they follow the period— K s magnitude relation indicated 
by the diagonal line in figure [2] Then, the amount of 
"deviation" from the K s magnitude that the period— K s 



magnitude relation predicts can be a measure of the de- 
gree of circumstellar extinction (Jal). For example, the 
observed K s band magnitude of an infrared star is more 
than 4 magnitude fainter than the prediction from the 
period— K s magnitude relation. The four magnitude ex- 
tinction in K s band correspond to A v of about 41 mag, and 
the optical depth at visual wavelength of about 38. Here 
we further assume that the extinction low for interstel- 
lar matter ( IU3I0 can be also applied to the circumstellar 
matter, although it is probably unrealistic. 

Under these assumptions, we can roughly estimate the 
mass loss rates of the infrared variables by using the es- 
timated A,, through the "deviations", van Loon ( fl2IO re- 
lated the total mass loss rate and optical extinction as, 
M[M e yr- 1 ] = 1.5 x lO- 9 Z[Z ]- a5 L[L ]°- 75 A°/ 75 . Sub- 
stituting typical LMC values of Z = 0.4, L = 8000, Av = 
20 mag gives M ~ 1.9 x 10~ 5 M Q yr-\ Av = 30 mag 
gives M ~ 2.6 x 10~ 5 M Q yr _1 , and Av = 40 mag gives 
M~3.2x 10~ 5 M @ yr- 1 . If we take L = 4000, instead 
of L = 8000, the mass loss rate is reduced by a factor of 
(4000/8000)° 75 - 0.6. In any case, it is likely that they 
are losing mass of the order of about 10~ 5 MQyr~ \ com- 
parable to the Galactic (i.e., presumably more metal-rich) 
infrared variables. These are only rough estimates, and 
the results of more detailed analysis will be published in 
the future paper. 

One more thing to note is that the A_gdeviationA.h 
from the period—^ magnitude relation and their ob- 
served near-infrared colors have quite tight relation Q@], 
and Ita et al. in preparation). This means that these near- 
infrared colors can be good indicators of mass loss rate 
for infrared stars. In other words, their degrees of circum- 
stellar extinctions can be estimated from their (J — K) or 
(H-K) colors. 
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